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Abstract 

We developed an easy, upscalable process to prepare lateral spin-valve devices on epitax- 
ially grown monolayer graphene on SiC(OOOl) and perform nonlocal spin transport measure- 
ments. We observe the longest spin relaxation times Ts in monolayer graphene, while the spin 
diffusion coefficient Ds is strongly reduced compared to typical results on exfoliated graphene. 
The increase of T5 is probably related to the changed substrate, while the cause for the small 
value of Ds remains an open question. 

Spin transport in graphene draws great attention since the observation of spin relaxation lengths 
of A5 = 2 jum, with spin relaxation times in the order of T5 = 150 ps at room temperature (RT) in 
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mechanical exfoliated single layer graphene (eSLG).'^ Recent experiments show an increase of T5 
to T5 ~ 0.5 ns in eSLG at RT^ ^ and T5 ^ 1 ns at T = 4 K.l^ Measurements on bilayer graphene 
(BLG) show even higher spin relaxation times, up to a few nanoseconds at low temperature.'^ At 
the same time, a study on few-layer graphene (FLO) showed an enhancement of T5 with increasing 
number of layers, which is attributed to the screening of external scattering potentials.!^ 
While most spin transport measurements were performed on exfoliated graphene, a first publi- 
cation by Avsar et showed measurements on graphene, grown by chemical vapor deposition 
(CVD) on copper foil. This publication marked the first step towards large scale production of spin 
transport devices, which showed similar spin transport properties compared to exfohated graphene. 
The disadvantage of the growth of graphene on metal substrates is however that one is forced to 
transfer the material to an insulating substrate to be able to perform transport measurements. 
Therefore it is useful to think about an alternative, e.g. epitaxially grown graphene on semi- 
insulating SiC.!^ This letter is the first report of spin transport in this material and therefore the 
first report of spin transport in graphene on a different substrate than Si02. We present lateral non- 
local spin-valve and spin-precession measurements on graphene strips prepared from monolayer 
epitaxial graphene (MLEG) grown on the Si-face of a semi-insulating SiC substrate (SiC(OOOl)) 
by sublimation of Si in Ar atmosphere.'^'^ 

4H-SiC wafers'^ are heated to 2000°C in an ambient argon pressure of 1 atm as described in Refs. 
[8] and IH leading to the growth of the so called buffer layer that is predominately (> 80%) covered 
with MLEG, with some areas uncovered or covered with double layer graphene. The measure- 
ments were performed on MLEG and with the help of Hall measurements on similar samples 
we estimate an electron doping with a charge carrier density of n ~ 3 x 10^^ cm^^ and a charge 
carrier mobility of /i ^ 1900 cm^/Vs at RT. 

Figure [Tj^a) shows an about 7x5 mm^ big part of a SiC wafer covered with MLEG, prepared 
with a pattern of Ti/Au structures that form a periodic pattern of bondpads with leads to central 
100 X 100 jUm^ areas for further device preparation. These structures are prepared in an optical 
lithography step, using a deep-UV mask aligner with a double resist layer (L0R-3A / ZEP-520A, 
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from MicroChem / ZEON Corporation). After development, the wafer is etched with oxygen 
plasma at 40 W for 20 seconds, before depositing a Ti/Au (5 nm/35 nm) double layer using 
e-beam evaporation followed by lift-off in PRS-3000 (from J.T. Baker). The etching step is nec- 
essary to enable the adhesion of the Ti/Au contacts on the substrate. To prepare the central device 
regions (Figure [TJ^b)), two MLEG strips per area are defined, using e-beam lithography (EBL) on 
a negative resist (ma-N 2400, from micro resist technology GmbH) and the uncovered MLEG is 
etched in a second oxygen plasma etching step. After this, the wafer is annealed for two hours 
in Ar(95%) : H2(5%) environment at 350°C to remove resist residues. To avoid the conductivity 
mismatch,'I^J^the wafer is covered with an approximately 1 nm thick AlOx layer by evaporating 
0.4 nm aluminum at a base pressure of ;? < 1 x 10^ mbar, oxidation in situ in O2 atmosphere at 
a base pressure of p > 3 x 10^^ mbar for 15 min and repeating the step a second time. Finally, in 
a standard PMMA resist based EBL step the 45 nm thick Co electrodes are formed (Figure [TJb)) 
connecting the Ti/Au leads with the two graphene strips, before the bondpads are contacted using 
wire bonding. By preparing the wafer with an optical step and using EBL only on the small central 
areas, we developed a fast and easy process to prepare a full wafer with (spin) transport devices. 
This process can also be used for different types of large area graphene on non-conducting sub- 
strates. 

As the resolution of the ma-N resist is limited to about 50 nm, we developed a second process, 
replacing the ma-N resist step with a PMMA step. This enables a higher resolution but in return 
requires an additional step to remove the graphene outside the 100 x 100 jim^ areas to disconnect 
the Ti/Au leads. This method was implemented for two devices, using an additional optical lithog- 
raphy step, where we cover the central device areas, while the exposed graphene is removed with 
oxygen plasma. 

The presented measurements are performed on a W = 0.7 /im wide MLEG strip in vacuum at 
a base pressure of about 1 x 10^ mbar using low frequency lock-in technique and AC currents 
between 1 and 10 /iA. The measurements have been confirmed with consistent results, that have 
been obtained on several spin-valve areas on two other devices with W = I /im. 
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Figure 1: (a) Optical microscope picture of a SiC wafer prepared with Ti/Au bondpads and leads 
to central 100 x 100 /im^ areas for further device preparation. Each of these patterns has a unique 
label for further production steps and measurements. The changes in the background color are due 
to scratches and resist residues on the backside of the transparent SiC wafer, (b) SEM picture of 
one of the central device areas with two spin valve devices, connected with Co electrodes to the 
Ti/Au leads, (c) Sketch of an MLEG spin valve device with four Co contacts. The wafer including 
the MLEG strip is covered with AlOx, before the Co contacts are deposited. 
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The typical nonlocal geometry is presented in Figure [TJc). A spin polarized current / is sent from 
contact 2 to contact 1, generating a spin accumulation at contact 2, that diffuses in positive and 
negative x-direction. The AlOx barrier separates the MLEG from the Co contacts and avoids re- 
absorption of the injected spins in the higher conducting cobalt.'^ The exponential decaying spin 
accumulation generates a nonlocal voltage Vni between the spin sensitive contacts 3 and 4, which 
can be measured as a function of the magnetic field. In a spin-valve measurement, the magnetic 
field By, aligned with the contacts, is first used to bring the magnetization of the electrodes into a 
parallel (P) configuration and is then ramped in the opposite direction. When the magnetization 
of one of the electrodes is switched the measured voltage shows abrupt changes. The magnetic 
switching fields of the contacts are different due to different coercive fields which are achieved by 
different width of the contacts .^^^ 

Figure|2];a) shows two spin-valve measurements, one at RT and one at 4.2 K. y„/, normalized to the 
nonlocal resistance = Vni/I, is plotted as a function of the magnetic field. The upper measure- 
ment has been obtained at RT. After saturating the magnetization of the contacts at 5^ ^ —450 mT 
no change of R„i is observed, before By crosses By = (not shown). Then at By = 18 mT a switch 
in Rni by 250 m^2 is observed, that can be attributed to the antiparallel alignment (AP) of injector 
and detector before it switches back to P and the initial value at By = 30 mT. Based on only 
two visible switches it can be concluded that the outer contacts (contact 1 and 4 in Figure ^c)) 
are giving no significant contribution to the signal.'^ The relatively small amplitude of the signal of 
2Rni ^ 0.3 t2 is not necessarily related to spin relaxation in the graphene strip but is here related 
to the relatively low polarization of the contact interface. Also, the contact interface is described 
by the R parameter, which is in our measurements R = W Rc/Rsq > 2.1 /im, with a contact re- 
sistance of 7?c > 3.3 kQ, a square resistance of the MLEG of Rsq = 1. 1 kf2 and W = 0.7 ^m^^^ 
Therefore the contacts are almost noninvasive, but can still slightly influence the spin transport 
measurements.'^^'^ 

The spin valve measurement, performed at T = 4.2 K, shows similar behavior. The main differ- 
ences are, that the amplitude is approximately doubled and the switching fields are slightly in- 
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creased, due to a change in the coercive fields with decreased temperature. Additional, shows a 
gradual decrease in its value, before the contacts switch to AP. This is probably due to a slight mis- 
alignment of the magnetic field and the electrodes. The changed background resistance is mainly 
influenced by heat related effects™ and can therefore change with temperature. 
To analyze the spin transport properties, we perform Hanle spin precession measurements.^^ For 
this purpose the magnetic field is aligned in z-direction. The resulting spin dynamics are described 
with the one dimensional Bloch equation for the spin accumulation fl^^ 



DsV^fis- — + (doxfis = d. 



(1) 



The first term on the left-hand side describes the spin-diffusion represented by the spin-diffusion 
coefficient Ds, and the second term describes the spin relaxation with the spin relaxation time T5. 
The third term describes the precession with the Larmor frequency cOq = gllB/h B, where ^ ^ 2 is 
the effective Lande factor and /ig is the Bohr magneton. 

The Hanle precession measurements in Figure[2];b)-(d) can be fitted with the solutions of the Bloch 
Table 1: Results of the fits to the measurements in Figure |2j;b)-(d). 



L(jum) 


2.9 


1.2 


1.2 


r(K) 


4.2 


4.2 


293 


D5(cmVs) 


4.06 ±0.05 


4.26 ±0.06 


2.38 ±0.03 


T5(ns) 


2.34±0.28 


1.66 ±0.02 


1.34 ±0.02 




0.98 ±0.06 


0.84±0.01 


0.56 ±0.01 




9.2±0.4 


275.7 ±2.2 


102.9 ±0.7 



equation ([T]), yielding the spin transport quantities Ds and T5. A summary of the fitting results are 
shown in Table [T] Figure [2];b) shows Hanle precession measurements, performed on a distance of 
L = 2.9 iim at 4.2 K. The fit gives Ds = 4.06 ±0.05 cm^/s and T5 = 2.34 ±0.28 ns, resulting in 
h = \fD^s = 0.98 ± 0.06 jum. 

We would like to note, that this value for T5 is the longest reported spin relaxation time on mono- 
layer graphene. And as the contacts are to some extent invasive, we can expect even higher values 
for T5, because a part of the injected spins relax at the contact interface. The effect of the contacts 
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Figure 2: Nonlocal spin transport measurements, (a) Spin valve measurements on a device with 
L = 1.2 jum inner contact distance at RT (purple and blue) and 4.2 K (red and green). The sweep 
directions of the magnetic field are indicated by the horizontal arrows, the relative orientation of the 
Co contacts is illustrated by the pairs of vertical arrows, (b)-(d) Hanle precession measurement for 
parallel (ft, blue boxes) and antiparallel alignment (t|, red open triangles) of the inner electrodes 
for (b) L = 2.9 /im at 4.2 K, (c) L = 1.2 ^um at RT and (d) L = 1.2 /im at 4.2 K. The fits to the 
solutions of the Bloch equation are plotted in gray. The background resistance, which is visible in 
(a), is subtracted in the Hanle precession measurements (b)-(d). 
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becomes apparent if one compares the fits of the measurement at L = 2.9 /im and L = 1.2 /im 
at 4.2 K (Figure [2]^b) and (d)). For the measurement at L = 1.2 /im we get T5 = 1.66 ±0.02 ns, 
which is around 70% of the T5 obtained from the L = 2.9 jUm measurement. This is due to the fact 
that the contact induced relaxation is more predominant, the shorter the distance the spins diffuse 
in the graphene strip between the contacts.— With invasive contacts, the shorter measurement dis- 
tance also leads to a slight increase in the measured Also this is observed, comparing the 
L = 1 .2 jUm to the L = 2.9 jUm precession at 4.2 K. 

When measuring at RT (Figure ^c)) we observe a reduction of Ds by more than 40% and T5 is 
decreased by about 20%. Therefore we get A5 = 0.56 ± 0.01 /im, which is one third smaller than 
Xs at 4.2 K (see Table [T]). We also observe slightly higher values for T5 of up to ~ 1.5 ns (not 
shown). Figure [3] shows a more detailed temperature dependence of T5, Ds, Xs and the non-local 
signal amplitude i?„/ between 4.2 K and RT. All four parameters show a decline between 4.2 K 
and RT. While Ds and A5 are monotonically decreasing by 40 and 30%, respectively, the value 
of Ts and R,ii drops by 20% and a factor of 3, respectively. The decrease of all four values can 
be related to electron-phonon scattering.'^^'^^' T5 and are approximately constant below 100 K 
which could be related to the fact, that the phonons are frozen out below that temperature. Given 
the relatively low mobility of the graphene, the temperature dependence could also be dominated 
by Coulomb scattering on trapped charges in the SiC substrate which shows a strong temperature 
dependence as described by Farmer et 

Our typical values for eSLG on Si02 at RT are in case of Ds about a factor of 80 bigger, however 
the measurements on MLEG strips show an increase of T5 by about a factor of 10. This still leads to 
a ~ 70% lower value for A^.'^^'^The increase in T5 in MLEG compared to eSLG can be attributed 
to the changed substrate. While Si02 has an electrical inhomogeneous surface potential leading 
to electron-hole puddles'^ and limiting effects for spin transport in graphene,^ the SiC crystal 
and the buffer layer are far more homogeneous and reduce therefore scattering. We would like to 
note that in our measurements on exfoliated graphene the spin transport properties are only weakly 
influenced by the temperature^ ^ whereas we here see an improvement at low temperatures. 
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Although Ts is improved, we do not know the origin for the reduced values of D5. We obtain 
the diffusion coefficient Ds from Hanle spin precession measurements. To verify if the value for 
Ds is correct, we compare it to the diffusion coefficient Dc acquired from charge transport mea- 
surements on the same area. Dc is calculated using the Einstein relation Dc = {Rsqe^v{Ef))^^ , 
where e is the electron charge and v is the density of states (DOS). The band structure for MLEG 

Ft! 

on SiC(OOOl) is the same as for eSLG. Therefore, we can assume the same DOS as for eSLG, 
v(£') = gvgs'^^ \E\ I {h\jr), with the twofold valley (^v = 2) and spin (gs = 2) degeneracies and 
the Fermi velocity Vf ~ 10^ ms^^. With n estimated by Hall measurements on similar devices 
and niEp) = jQ'v{E)dE we can calculate the Fermi energy Ef and receive v{Ef)- With n ^ 
3 X 10'^ cm^^ and Rsq = 1.1 kQ. we get Dc ~ 190 cm^/s which is similar to typical values ob- 
tained in eSLCl^This is not surprising, because the charge carrier mobility in our samples is with 
/i = (Rsqen)^^ ^ 1900 cm^/Vs reasonably close to the mobility of our eSLG devices.!^ 
But this value of Dc means that we observe a difference between the charge diffusion coefficient 
Dc and Ds of a factor 45 to 80 (compare Table[T]) in contrast to Dc ~ Ds in eSLG. and Dc can 
in principle be different as observed by Weber et al.l^^lin a two-dimensional electron gas. Here the 
effect is attributed to electron-electron interactions but was significantly smaller than in our results. 
In FLG a slight difference of ~ 20% between the two coefficients is found, ^ but that difference is 
not comparable to the observation here. 

We therefore do not expect a difference between the diffusion coefficient obtained from charge 

transport measurements and from Hanle precession measurements. While we cannot explain the 

observed difference, yet, we can exclude some possible explanations for it. 

We do not expect the Dc value to be incorrect as the observed charge transport is comparable to 

eSLG. One aspect though, that could result in a wrong Dc value, are extra current paths next to the 

MLEG strip which would result in a change of the observed R^q. We can exclude this by carefully 

controlling the etched structures with a scanning electron microscope (SEM) and by confirming 

that contacts of different MLEG strips show no conduction between each other. 

n was determined by Hall measurements on similar samples, but not on the spin transport samples 
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Figure 3: Temperature dependence of (a) the spin relaxation time, (b) the spin diffusion coefficient, 
(c) the spin relaxation length and (d) the nonlocal signal for the sample with L = 1.2 jum. If 
available, several fitting results at the same temperature were averaged. 
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themselves, therefore there could be an error in the value for n which would lead to an incorrect 
value for the DOS. The highest values for n, measured on MLEG samples on SiC(OOOl), under the 
described growth conditions, are around I x 10^^ cm^^ which leads to Dc ~ 100 cm^/s. This 
changes Dc by less than a factor of 2 and smaller values for n would only increase the calculated 
value for Dc- Hence, also this aspect does not explain the difference between Ds and Dc- 
Another possibility would be a wrongly assumed DOS. Though, to result in values for Dc similar 
to Ds, we would need a DOS as high as ~ 50 times the DOS of BLG. But we can be sure that we 
do not have such a DOS in our system because similar material to that used in our studies shows 
the typical quantum Hall effect (QHE) of eSLG.ISES] 

Since we do not find any explanation for the difference between Dc and Ds in the way Dc is deter- 
mined, let us have a look at Ds- Ds is obtained by the fit of the Hanle precession data in the same 
way as in earlier experiments J ^ 1 ^ I ^ ^ I S I Therefore the fitting procedure cannot explain the difference 
in the values, as the result for eSLG, Dc ~ Ds, is based on fits to measurements on eSLG using the 
same procedure. Next to that, we can confirm the value for Ds in a different way. As mentioned 
before, the small value for Ds leads to a relatively small value for A^. The order of magnitude 
of this value can be confirmed by comparing the change of with L. By fitting an exponential 
decayf^to the two values versus L for the data obtained at 4.2 K (see Table [T| fit not shown), 
we receive A5 ~ 0.5 jim in agreement with the order of magnitude of A5 obtained from our fitted 
Ts and Ds^ With Ds = Dc ~ 200 cm^/s and T5 ~ 2 ns we would receive a A5 of one order of 
magnitude larger. 

While Xs is confirmed, there is still the possibility that the prefactor of the precession term in the 
Bloch equation ([T]) is wrong, which would affect linearly the determination of Ds and 1 / T5. This 
would be the case, if the effective Lande factor g is changed in our system. But the reduction of g 
by a factor of ~ 50 is needed, to result in our measured Ds- This is unlikely, also since ^ ~ 2 was 
confirmed for epitaxial multilayer graphene on the C-face of SiC^ 

A change of Ds can only be caused by the substrate as we expect the graphene to be comparable 
to cSLG'^ and growth related defects like grain boundaries do not show a strong effect on Ds and 
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Ts for CVD grown single layer graphene on 8102- One of the substrate related effects could be 
inhomogeneities of the graphene thickness and doping at terrace step edges'^ and scattering poten- 
tials resulting from that. However, this is relatively unlikely since step edges are not resulting in a 
discontinuity of the graphene layer.l^l On the other hand, the out of plane electric field between the 
bulk SiC and the buffer layerl^ could have an effect on the spin transport. 

Another possible reason for the change in the spin transport properties could be related to the 
buffer layer. Its topology is graphene-like, though a part of the C atoms is covalently bonded to 
the underlying Si atoms. Therefore, the layer is electrically inactive and only weakly interacts 
with graphene layers on top.l^ The buffer layer does not seem to affect charge transport, at least 
not the measured resistance or the QHE,^ although it influences the temperature dependence of 
the charge carrier mobility.!^ However, localized states in the buffer layer could act as hopping 
sites for electron spins and could influence the spin relaxation and the spin diffusion. By spins 
hopping into these states and back, Ds could be reduced without affecting Rsq and therefore the 
determined Dc (as we do not include these extra states in the DOS). This kind of localized states 
could also originate from Al clusters in the AlOx barrier. When depositing the barrier, some of 
the Al atoms could cluster and if their size exceeds some certain limit, part of the Al could stay 
non-oxidized. Those clusters would have a high DOS compared to the MLEG and could there- 
fore have a relatively strong influence on the diffusion. This effect is unlikely, as we do not see it 
for eSLG on Si02 but the less rough surface of MLEG on SiC and the resulting difference in the 
growth mechanism could result in this clustering. Here it would be interesting to produce samples 
with the AlOx barrier only locally below the contacts as discussed in Ref. [14] to compare the spin 
transport properties with the here reported results. 

One other effect that could affect the measurements is the influence the Ar(95%) : H2(5%) clean- 
ing at 350°C could have on the buffer layer. F. Speck et al. discuss the intercalation of hydrogen 
in epitaxial graphene on SiC which leads to the transformation of the buffer layer into an extra 
graphene layer. Though the discussed experiment uses about 1 bar pure hydrogen at 550°C for 
75 minutes, our cleaning step could partly intercalate hydrogen below the graphene layer and this 
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could lead to extra transport channels and influence the transport measurements. 
Non of these considerations above led to a conclusive explanation of the observed difference be- 
tween the diffusion coefficients obtained from charge and spin transport measurements. Therefore 
we have to wait for further measurements to determine if the difference is based on the way those 
values are obtained or if there is a difference between charge and spin diffusion in MLEG on 
SiC(OOOl). The effect of the buffer layer can be addressed by measuring spin transport on quasi 
freestanding MLEG on SiCp^and general substrate related effects can be examined by transferring 
MLEG to SiOi^ 

In summary we present a fast and easy process to prepare (spin) transport devices on wafer scale 
graphene by the example of MLEG. With this technique we produced lateral spin- valve devices on 
MLEG and performed spin-valve and Hanle spin precession measurements between T = 4.2 K and 
RT. In the Hanle measurements we observe exceptionally high values for T5 of up to T5 = 2.3 ns 
and very small values for Ds of Ds < 5 cm^/ s, resulting in a reduction of A5 by a factor of 2 to 3 
compared to eSLG. We observe a significant difference between the diffusion coefficient obtained 
from charge and spin transport measurements, which we discuss but cannot explain, yet. Finally 
we present the temperature dependence of the spin transport and show a decrease for T5, Ds, A5 
and Rfii with increasing temperature, that can be linked to electron-phonon scattering or Coulomb 
scattering on trapped charges in the SiC substrate. 
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